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This invention relates to a signal processing apparatus and a method of 
signal processing. More particularly, but not exclusively, the invention 
5 relates to a signal processing apparatus and method of signal processing 
for use in multiple input multiple output (MIMO) antenna architecture. 

Current iterative processing techniques allow wireless telecommunications 
data transmission rates for. links between single antenna transceivers that 
10 approach the Shannon data capacity limit, for a single input single output 
(SISO) system. The Shannon data capacity limit is given by: 



Where: 

15 C is the channel capacity 

W is the channel bandwidth in Hz 
S/N is the signal-to-noise power ratio. 

The concept of spatial multiplexing allows the realisation of enormous 
20 data transmission rates when multiple transmit and receive antennas are 
employed. For a given average signal-to-noise ratio (SNR) per receive 
antenna, or total transmit power, the theoretical MIMO capacity increases 
linearly with the number of transmit or receive antennas, whichever is the 
smaller, at high SNR values. Therefore, it is desirable to implement 
25 systems with large MIMO configurations, for example four or more 
transmit antennas (Tx) and four, or more receive antennas (Rx), or more. 
These large MIMO configurations are typically integrated with iterative 
processing in order to provide reliable transmission of data at high 

■ 

spectral efficiencies, an example of such a prior art MIMO architecture 
30 combined with iterative decoding is shown in Figure 1. 



C = W log 2 (1 + S/N) bit/s 
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One problem associated with large MIMO architectures is a high 
computational burden that is imposed upon them by the execution of Log 
Likelihood Ratio (LLR) calculations at a front end of an iterative decoder 
on the receive side of the MIMO architecture. This is due to the LLR 

5 calculations having to account for every possible candidate symbol across 
every transmit antenna. These calculations require marginal probabilities 
to be obtained by summing over all symbol probabilities to obtain bit LLR 
values. The number of these calculations grows exponentially with the 
number of transmit antennas and the modulation level. The number of 

10 calculations that must be executed to obtain LLR value could be 
thousands, tens of thousands or even up to millions per second. 

According to a first aspect of the present invention there is provided a 
signal processing apparatus comprising a plurality of receiving means 
15 arranged to receive a composite signal indicative of a plurality of symbols 
transmitted, substantially simultaneously, from a plurality of remote 
transmission means, and processing means arranged to iteratively decode 
each most probable symbol contained in said composite signal, within a 
constrained enumeration formalism. 

20 

The processing means may be arranged to define an enumeration 
constraint for use in the constrained enumeration formalism. 

The use of an enumeration constraint in such an apparatus reduces the 
25 computational load required to determine the bit stream when compared to 
prior art arrangements. This either frees up processing power for other 
applications or allows a reduced computational and/or electrical power, 
cheaper processor to be used in determining this bit stream. The 
reduction in the computational load also increases battery life in portable 
30 equipment, as the processor need not be operated as frequently in prior 
art arrangements. 
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The processing means may be arranged to perform a QR decomposition 
upon a channel gain matrix. The enumeration constraint may be a number 
of entries in an R matrix over which probable symbols are enumerated. 
5 This establishes a small area within an R matrix for symbol enumeration 
rather than using a whole channel gain matrix, thus limiting the size of a 
list of possible candidate symbols. The enumeration constraint may be of 
constant dimension. 

10 The processing means may be arranged to determine the most probable 
symbol by enumerating across all symbol conditional probabilities for 
each possible symbol. 

The processing means may be arranged to convert symbol conditional 
15 probability to a bit level, or a symbol level, logarithmic likelihood ratio 
(LLR) . 

The apparatus may include a parallel to serial conversion means arranged 
to convert parallel, bit level, LLR's into a single stream of LLR's. 

20 

The apparatus may include a deinterleaving means arranged to 
deinterleave, bit level, LLR's from the single stream of LLR's. This 
improves data channel performance as it improves the assumption of 
independent transmitted bits. 

25 

The apparatus may include decoding means arranged to apply iterative 
soft input soft output (SISO) decoding single bit LLR's to determine a 
symbol. The SISO decoding may be a Maximum-A-Posteriori (MAP) 
technique or it may be a Low Density Parity Check (LDPC) technique. 



30 
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The decoding means may be arranged to pass a symbol, or bit, probability 
to the processor for inclusion in an iterative enumeration step. 

The apparatus may include a hard decision unit that is. arranged to 
determine a symbol based upon a soft output from the decoding means. 

According to a second aspect of the present invention there is provided a 
method of signal processing for a MIMO system comprising the steps of: 

i) receiving a composite signal indicative of a plurality of symbols; 

ii) performing a QR decomposition upon a channel gain matrix for 
the composite signal; 

iii) defining an enumeration constraint; 

iv) calculating possible conditional probabilities for one of the 
plurality of symbols contained within the composite signal, using 
the enumeration constraint; and 

v) iterating step iv), incorporating a most probable symbol for the one 
symbol determined in the previous iteration of step iv) in the 
conditional probability calculation operation. 

The method may include setting the enumeration constraint to encompass 
a sub-set of possible transmit antennas. 

The method may include defining the enumeration constraint as a number 
of elements within an R matrix, for example an enumeration sub-matrix 
within the R matrix. 

The method may include calculating a symbol conditional probability in 
order to determine the most probable symbol transmitted by a given 
transmit antenna. The method may include converting a symbol 
conditional probability to a bit level logarithm likelihood ratio (LLR). 
The method may include calculating a symbol conditional probability in 
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order to determine the most probable symbol received over a given 
transmission channel. 

The method may include converting a plurality of parallel streams of bit 
5 level LLR's to a serial stream of bit level LLR's. 

The method may include deinterleaving bit level LLR's from the serial 
stream of bit level LLR's. 

10 The method may include decoding the single bit LLR's, typically using a 
soft input soft output (SISO) technique, such as for example, a Maximum- 
A-Posteriori (MAP) technique or a Low Density Parity Check (LDPC) 
technique. 

15 The method may include making a hard determination of a received 
symbol based upon a soft output from the decoding step. 



According to a third aspect of the present invention there is provided a 
method of reducing the computational load of a signal processor in MIMO 
20 architectures comprising the steps of: 

i) receiving composite input signals having spatial diversity from 
each of a set of n receiver elements; 

ii) constructing an n by m channel matrix from values indicative of 
channel gains between each transmit and receive element; 

25 iii) executing a QR decomposition upon the channel matrix to form an 

upper triangular R matrix and a unitary Q matrix; 

* 

iv) enumerating to determine probabilities of a given symbol being 
transmitted from a given transmitter using a constrained data sub- 
set of the triangular matrix; and 
30 v) making a hard decision about which possible symbol is the most 

probable symbol to have been transmitted so as to reduce the 
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number of enumerations required to carry out a further probability 
calculation. 

The method may include using sub-optimally determined symbol values to 
5 generate final definite symbol values. 

« 

The method may include using the most probable determined symbol to 
generate a final hard symbol value. 

10 According to a fourth aspect of the present invention there is provided a 
computer readable medium having therein instructions for causing a 
processing unit to execute the method of either of the second or third 
aspects of the present invention. 

15 The invention will now be described, by way of example only, with 
reference to the accompanying drawings, in which: 

Figure 1 is a schematic representation of a MIMO architecture of 
the prior art; 

20 

Figure la is an idealised representation of a constellation diagram 
for 16 QAM; 

Figure lb is a representation of a the constellation diagram for 16 
25 QAM between a transmitter and receiver with a complex channel 

gain; 

Figure 2 is a schematic representation of an embodiment of a 
MIMO architecture in accordance with an aspect of the present 
30 invention; 
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Figure 3 is a schematic representation of components a joint 
detect-decode unit of Figure 2; 

Figure 4a is a schematic representation of an initial step of 
5 enumeration across possible symbols sent from transmit antennas 3 

and 4 of a four transmit antenna, four receive antenna MIMO 
architecture; 

Figure 4b is a schematic representation of a step of enumeration, 
10 subsequent to that of Figure 4a, across possible symbols sent from 

transmit antennas 2 and 3 of the four transmit antenna, four receive 
antenna MIMO architecture; 

■ 

Figure 4c is a schematic representation of a final step of 
15 enumeration across possible symbols sent from transmit antennas 1 

and 2 of the four transmit antenna, four receive antenna MIMO 
architecture; and 

Figure 5 is a flow diagram detailing a method of signal processing 
20 in accordance with an aspect of the present invention. 

Referring now to Figure 1, a typical MIMO architecture 100 comprises a 
transmission device 102 and a reception device 104. 

25 The transmission device 102 comprises an encoder 105, an 
interleaver 106, a serial to parallel converter 108 and mappers 110 a-d, 
each mapper 110 a-d having a corresponding transmit antenna 112 a-d. 

The encoder 105 channel codes source data using, for example Parallel 
30 Concatenated Convolutional Codes (PCCC), Serial Convolutional Codes 
(SCC) or Low Density Parity Checks (LDPC) . The channel coded source 
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data is then interleaved by the interleaver and divided into m parallel, 
independent data steams at the serial to parallel converter 108, in this 
example four data steams. ' The four parallel, independent data streams 
are converted to complex symbols having a plurality of bits, typically, 
5 according to Gray coded constellation mapping, by the mappers 110 a-d. 
Thus, this transmission strategy incorporates Bit Interleaved Coded 
Modulation (BICM). An individual complex symbol is transmitted, 
simultaneously, by each of the transmit antennas 112 a-d. 

10 The transmitted signals pass across independent complex Gaussian 
distributed channels, typically in a rich scattering environment containing 
multiple, randomly arranged scatterers that may be stationary or mobile, 
to the reception device 104. 

15 The reception device 104 comprises n,' in this instance four but it could be 
more, or less, than four, receive antennas 114 a-d, channel state 
information (CSI) estimators 116 a-d and a demapper 118. The reception 
device also includes a parallel to serial converter 120, a deinterleaver 122 
and a decoder 123. 

20 

Each receive antenna receives all of the symbols simultaneously 
transmitted from the transmit antennas via a MIMO channel. Thus each . 
of the receive antennas 114 a-d receives a composite signal containing 
superposed weighted components of each of the transmitted symbols. The 
25 specific weighting of each of the symbols is given by a channel gain value 
for the channel between a given transmit antennas and a given receive 
antenna. 

If the individual channels between each given transmit antenna and each 
30 receive antenna is different each receive antenna will receive a different 
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composite signal. It is this spatial diversity upon receive that allows a 
symbol determination to be achieved. 

* 

In order to determine what symbols were transmitted across the channels 
5 prior art systems determine the probability of a given symbol being 
transmitted from a given transmit antenna by marginalising the other 
transmitted symbols. 

■ 

Consider the case of a four transmit antenna, for four receive antenna 
10 MIMO system transmitting 16-QAM (Quadrature Amplitude Modulation) 
symbols, i.e. symbols from a set of 16 possible symbols. In order to 
marginalise three of the transmitted symbols it is necessary to cycle 
through the set of sixteen possible symbol values for each of the three 
transmitted symbols, as shown below. This requires 16 3 probability 
15 calculations to be executed in order to determine the marginal probability 
of a transmitted symbol, C jt being a specific value, S p . This process 
must be repeated for the given symbol having each of the sixteen possible 
symbol values leading to a total computational load of 16 4 probability 
calculations. 

20 

P iW = I p[y>y 2 y>y<\c x c 2 c,c a ] 

S p z{S x ,S 2 , ,S l6 } 

PtCj^S!) is only concerned with probability of the symbol, S,, having 
been transmitted from transmit antenna 1, it is not concerned the symbols 
transmitted from the other three antennas of the systems. 

25 

In order to calculate the probability of symbol S, having been received it 
is necessary to marginalise out the probabilities of the contributions from 
the other three transmitted symbols. This is achieved by summing the 
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probabilites of all possible combinations of the other three transmitted 

■ 

symbols whilst holding Q = Si: 
P[Y|, Y 2 , Y3, Y 4 |Si, Si, Si, Si] 



5 P[Y„ Y 2 . Y 3> Y«|S lf Si, S„ S 16 ] 
P[Yi, Y 2 , Y 3 , Y 4 |Si, Si, S 2 , Si] 



P[Y if Y 2 , Y 3 , Y< 1 Si, Si 6 , Si 6f Si 6 ] 



10 This must now be repeated to calculate the marginal probabilities for 
d = S 2 to S l6 in order that all symbol probabilities are calculated. The 
SISO operates at bit level LLR values and therefore conversion from 
symbol probabilities to bit level LLR values is required. 



15 An example of a LLR calculation for a single bit is: 



LLR (6,)= log 



Prft, =1) 
Pr(6,=0) 



for a single bit 



Where Pr (bj = 1) is the probability of a given bit, b, having a value 1 
20 Pr (b,= 0) is the probability of a given bit, b, having a value 0 

Thus when Pr (b, = 1) = 1 and Pr (b, = 0) = 0 LLR (b,) -> + co 
Similarly when Pr (b, = 1)= 0 and Pr (b, = 0) = 1 LLR (b,) -> - oo 

25 And when Pr (b, = 1) = Pr (b, = 0) = Vfc LLR (b,) 0 



WO 2004/088908 PCT/GB2004/001386 

t 

11 

In the general case for symbols, the probability of a symbol transmitted 
from a given antenna, y'\ being one of a set of transmitted symbols, C at 
the given transmit antenna, y, where 1 <j<m and S p e {Sj,S2,: ,S 2 m c } - 

The sum of conditional probabilities, conditioned on the transmitted 
5 symbols, summed over all possible transmitted symbol combinations from 
all other transmit antennas is given by: 

= s p )= £/>[y|c] 

10 The whole process must be repeated for C 2) C 3 and C 4 thus leading to a 
- very computationally intensive process for the determination of received 
symbols. 

Once the codewords, C u C 2 ,C 3 , C 4l have been determined they are 
15 compared to codewords in a codebook, typically a look up table, in order 
to determine what bit sequence each of them correspond to, and this bit 
sequence is output from the receive unit. 

However, considering the case of 16QAM, and idealised received 
20 codeword constellation diagram consisting of all possible 16 symbols (SI- 
S16) in 16 QAM is show in Figure la. In the case of a hard decision 
making system where symbol S7 was transmitted, systems noise, for 
example Johnson noise or shot noise within the receiver, can cause the 
received signal to deviate from the constellation ordinate corresponding to 
25 symbol S7 towards, for example, the ordinate corresponding to symbol 
S3. When such deviation is small a correct hard decision-, that the 
transmitted symbol is symbol S7 is made. However, when such deviation 
is sufficiently large an incorrect hard decision is made that the transmitted 
symbol is symbol S3 rather than symbol S7. This is because there is a 
30 higher probability that the transmitted symbol is symbol S3 rather than 
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symbol S7 in this case as shown by the received symbol lying closer to 
symbol S3 than symbol S7 in the constellation diagrams. Such incorrect 
determinations of a symbol clearly leads to an error in the bit streams 
output from the receive unit. 

5 

In the case of a soft decision systems an incorrect codeword may still be 
determined based upon the basis of a higher calculated marginal 
probability however, the determined codewords are compared to 
allowable sequences of codewords prior to conversion to output bits. 
10 This comparison allows non-allowable codeword sequences to be 
reviewed and errors to be corrected before conversions to an output bit 
level data stream. 

In the physical case, each receive antenna will have an independent 
15 distorted symbol constellation for symbols received from a given transmit 
antenna as the channel between transmit and receive antennas is 
independent, see for example Figure lb. The effect of additive, white 
Gaussian noise, for example shot noise and Johnson noise, is to cause 
time varying random variations in the location of received signals (Si-Si 6 ) 
20 away from the ideal symbol location within a constellation. A function of 
Gaussian noise is that the mean of the distribution is the mean, thus the 
mean locations of the time varying received signal within the constellation 
is the actual locations of the transmitted symbol. A further feature of this 
noise is that it is symmetric about its mean point and is also independent 
25 of the receive antennas. 

> 

As the statistical processes, i.e. the additive Gaussian white noise, that 
cause the deviation of the received symbol away from its ideal 
constellation position are independent of each other the probabilities of a 
30 given signal being transmitted from a given .transmit antenna, and ■ 
received at the receive antennas can be multiplied in order to yield a joint 
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probability of a given set of signals being transmitted and received across 
all receive antennas, as follows: 

' C:Cj-S„ 

S,efo,S 2 , S i6 ] 

5 

■ 

Using current techniques it is still necessary to enumerate through all 
possible 16 4 values of C for each receiver as any variation in C causes a 
change in the mean received symbols at each receive antenna. 

10 Referring now to Figures 2 and 3, a MIMO architecture 200, of at least 
an aspect of the present invention, comprises a transmission device 202 
and a reception device 204. The transmission device 202 is substantially 
similar to that described in relation to Figure 1 and similar parts are 
accorded similar reference numerals in the two hundred series. 

15 

Typically, signals transmitted from the transmission device 202 are 
transmitted across a scatter rich transmission environment 215, which 
introduces spatial diversity upon receive into the signal received at the 
reception device 204. 

20 

4 

The reception device 204 comprises a four of receive antennas 214 a-d, 
channel state information (CSI) estimators 216a-d, a combined joint 
detector-decoder 218 and a hard decision making device 220, a slicer. 

■ * 

25 The receive antennas 214 a-d are arranged to receive the signals 
transmitted from the transmit antennas 212 a-d. The signals are passed to 
the joint detector-decoder unit 218 where statistics, usually in the form of 
a log likelihood ratio (LLR), are calculated to produce a likelihood of a 



• 
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particular symbol being present at a given point in a data stream, as will 
be described hereinafter. The statistics are passed to the hard decision 
making device 220 where they are processed to make a hard decision as to 
which symbols are present within the data stream from a set of possible 
5 transmitted symbols; 



The joint detector-decoder unit 218 comprises a QR decomposition 
unit 221, which receives inputs from the receive antennas 214 a-d and 

■ • 

■ 

also channel state information (CSI) 222 from the CSI estimators 216 a-d, 
10 a symbol probability processor 224, a symbol to bit LLR converter 225, a 
parallel to serial converter 226, a deinterleaver 228 and a SISO 
decoder 230. 



The operation of. the QR decomposition unit 221 is described immediately 
15 hereinafter. 



A MIMO data channel can be represented by: 



20 



Y = HC + 77 



Where H = 



H \\ H \2 H iZ H IA 

# // 2l H 22 H 22 H 2A 



H 



44. 



for a four transmit antenna, four receive antenna system. 



25 



H U H\l H \Z H \A 

H 2] H 22 H u H u 
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where: 

Y h Y 2 , Y 3 and Y 4 are the received signals in a composite signal received at 
the four receive antennas. 

5 For the general case K, is the receiving signal at receive antenna i where 1 
<i <n and n is the number of receive antennas. 

C is constructed from a matrix of mM c bits, where m is the number of 
transmit antennas and M c is the number of bits per symbol. For example, 
10 for 16 quadrature amplitude modulation (16-QAM) M c = 4. 

■ 

The matrix C for the four transmitter 16-QAM case is given by C, 



b = 



b n 


b \l 


b l2 


bi 4 


b 2 \ 


b ll 


b 23 


b 14 


6 31 


hi 


b 33 


634 


b A\ 


b 42 


643 


&44 

•< 



15 C = map (b) = [C, C 2 C„] 

C, to C M are transmitted symbols 

77 is a column vector of complex additive white Gaussian noise (AWGN) 
samples taken from a Gaussian distribution with zero means and variance 
20 per dimension, a 2 and is a function that is representative of noise within 
the system, both general thermal noise and component noise. 

The channel gain matrix, H, can be described by an orthogonal matrix, 
Q, and a matrix, R, having an upper triangular structure, subject to the 
25 channel gain matrix being a full rank matrix, therefore, there is sufficient 
diversity within the channels. Such a description of the channel gain 
matrix is called a QR decomposition of the channel gain matrix. 
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#11 #12 #13 # 



14 



#21 #22 #23 #24 



H 



44 _ 
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Qn Qn Qn Q 



■ Q 



44 _ 



K.1 


R l2 


*,3 




0 




"^23 


^24 


0 


0 


^33 


*34 


0 


0 


* 

0 


^44. 



As Q is orthogonal Q T Q = I as Q T = Q 1 
Q >.H = R 

Y = HC + r\ 
= QRC + Ti 
therefore Q'Y = RC + Q't] 

Both Q' 1 and R are fixed, and R has fewer non-zero elements when 
compared to H as the lower left hand side of R is zeros. 

Q'Y = Y = RC + fj 



It should be noted that the statistics of the noise vector, T|, are invariant 
under the transformation of Q" 1 due to their being circularly symmetric. 



In the general case and accounting for any previous decision made 

■ 

15 relating to transmitted symbol values the marginalised, conditional 
symbol probabilities are given by: 



20 



f(Cj = C p ) = 



A 

Ci 



2a' 



where E c <x <m with the constraint E c limiting the vector C to C 

[^2 + ] " Ec» C*-i2 " Ect iC x ] 



m 



The third term of the numerator , C, , is not included in an initial 



decision for C m , based upon {C iip ) as this term describes interference 
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cancellation based upon previously made symbol decisions. 

■ 

As Y is distributed by fj in the same way as Y is by *n. Y can be 
described as follows: 



5 



Y\ 




"*„ 




*,3 








Y 2 




0 


R 22 


R 2J 


R 24 






Yi 




0 


0 


^33 


^34 






74. 




_ 0 


0 


0 


^44. 




C4. 



As can be seen from the above 7 4 is dependent only upon the matrix 
element R 44 . Thus, using the above arrangement it is possible to 
determine the most probable symbol, S p , by enumerating only through, in 
10 the case of 16QAM, the sixteen possibilities for the value of C 4 . Whilst 
this is theoretically possible it would only be acceptable when the signal 
to noise ratio is very high. This is because if the mean value of a specific 
channel gain, in this case represented by R 4 „, is small the effect of noise 
in disrupting the channel will be significant. 

15 

It is usual to enumerate through the possible symbols over two channels, 

for example Y y and Y 4 in order to ensure that there is sufficient diversity 

within the enumerated channels so as to reduce the possibility of all three 
channel gains enumerated being small enough to be affected by noise. 
20 This arrangement reduces the complexity of determining transmitted 
symbols as, in the example discussed above, d and C 2 are not required to 
determine C 3 and C 4 . 

Referring now to Figure 4a, in the example shown below there are 
25 assumed to be four transmit antennas, m, and four receive antennas, n, 
and an enumeration constraint 401, E c , of two, i.e. only two symbols are 
considered for determination at any one time. 
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The conditional symbol probability including Gaussian noise, is given by: 



P(C 4 = S D )cc 



1 




z ri ex p 



271C 2 I " C:C 4 =S- »=3 



*=3 



2o 



when this calculation is dependent solely upon transmit and receive 
antennas 3 and 4. 



10 Now C is limited to C = [C 3 ,C 4 ], as shown in Figure 4a, the upper 

right-hand side of the matrix R contains the channel information values, 
represented as squares 402 in Figures 4a and the lower left-hand side 
represents null, zero 404, values of the matrix. 



15 The symbol with the highest conditional probability is selected by the 
symbol probability processor 224 as an estimated transmitted symbol 
transmitted from transmit antenna 4 and is given by: 

In effect the joint probability calculation is simplified to : 

Pte^r S *|c]p[nfe] 

C:Cj~S„ 

^e{5„S 2 , 5 I6 } 

20 



by the exclusion of C, and C 2 from consideration. 
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Once C, is determined as S n this value can be fixed for all subsequent 
calculations as this merely affects the mean value of subsequent 
calculations: 



X 






^12 


^13 








0 


^22 




^24 


y, 




0 


0 




R» 


A 




0 

-• 


0 


0 


^44. 



5 In an alternative computational scheme the fourth column of R is 
multiplied by S 0 and subtracted from Y thereby removing the effect of the 
previously determined from Y . 
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The effect of fixing C 4 = S„ is shown in Figure 4b wherein the marginal, 
conditional probability of any given symbol being transmitted from 
transmitter 3 is then calculated according to: 



P(C 3 =S p )«: £ fjexp 



C:C 3 =S p i = 2 



r.-t 



k = 2 



2a 



As can be seen this conditional symbol probability calculation takes into 

A 

account the previous symbol estimate made for Ca and transmit and 
receive antennas 2 and 3 are considered for this calculation. The third 
term in the numerator cancels out the interference from transmit 
antenna 4, shown in Figure 4b as a solid square 406. 



A soft decision is now made for the symbol transmitted from transmit 
20 antenna 3 based upon the conditional probability of the most likely 

■ 

A 

symbol to have been transmitted, C, where j = 3, in a similar manner as 
that for the symbol transmitted from antenna 4 as described hereinbefore. 
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Referring now to Figure 4c conditional probabilities of the symbols 
transmitted by transmit antennas 1 and 2 to receive antennas 1 and 2 are 



now calculated 



P(Cj=S p )* ]T JJexp 



C:Cj=S p i=l 



2 4 A 

*=1 1 = 3 



\ 



2a 



The vector C = [C!,C 2 ] is used to represent these conditional probabilities. 
The third term in the numerator cancels out the interference effects from 
antennas 3 and 4, shown in Figure 4c as solid squares 406. 



10 Again the transmitted symbols are estimated as being the symbol with the 
largest calculated marginal, conditional probability. 



The estimated symbols probabilities are converted to conditional bit 
probabilities at the symbol to bit LLR convertor 225 through the 
15 following: 



Cj:Cj=m»p(b / , / -{!,.0]) 

where for 16-QAM: 
20 bjb 2 b|bo 
P(C t )= P(0000) 
P(C 2 ) = P(0001) 
P(C 3 ) = P(0010) 
P(C<) = P(0011) 
25 P(C 5 ) = P(0100) 
P(C 6 ) = P(0101) 
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10 



P(C 7 ) 
P(C 8 ) 
P(C 9 ) 

P(C I0 ) 

P(C„) 

P(C 12 ) 

P(C 13 ) 
P(C 14 ) 

P(C 15 ) 
P(C 16 ) 



P(0110) 
P(0111) 
P(1000) 
P(lOOl) 
P(1010) 
P(1011) 
P(llOO) 
P(1101) 
P(1110) 

P(llll) 



Thus 



15 P(b 0 = 


0) 


= P(C.) 


+ P(C 3 ) + P(C 5 )+P(C 7 ) 


+ P(C 9 ) + P(C M ) + P(C 13 ) + P(C l5 ) 


P(b 0 = 


1) 


= P(C 2 ) 


+ P(C 4 ) + P(C 6 )+P(C 8 ) 


+ P(C l0 ) + P(C 12 ) + P(C 14 ) + P(C 16 ) 


P(b,= 


0) 


= P(C.) 


+ P(C 2 ) + P(C 5 ) + P(C 6 ) 


+ P(C 9 ) + P(C,o) + P(C 13 ) + P(C 14 ) 


P(b,= 


1) 


= P(C } ) 


+ P(C 4 ) + P(C 7 ) + P(C 8 ) 


+ P(C ll )+P(C 12 ) + P(C l5 ) + P(C l8 ) 


P(b 2 = 


0) 


= P(C.) 

• 


+ P(C 2 ) + P(C 3 )+P(C 4 ) 


+ P(C 9 )+P(C 10 ) + P(C 11 ) + P(C 12 ) 


20 P(b 2 = 

• 


1) 


= P(C 5 ) 


+ P(C 6 ) + P(C 7 ) + P(C 8 ) 


+ P(C 13 )+P(C 14 ) + P(C 15 ) + P(C l6 ) 


P(b 3 = 


0) 


= P(C.) 


+ P(C 2 ) + P(C 3 ) + P(C 4 ) 


+ P(C 5 ) + P(C 6 ) + P(C 7 )+P(C 8 ) 


P(b 3 = 


1) 


= P(C 9 ) 


+ P(C,») + P(C„) + P(C 12 ) + P(C l3 ) + P(C l4 ) + P(C l5 ) + P(C 16 ) 



25 



Thus, bit level log likelihood ratios (LLR's) can be given by the natural 
logarithm: 



A(6 y ,/) = log 
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The bit level LLR's are fed into the parallel serial converter 226 where 
they are converted into a single stream and are then passed to the 
deinterleaver 228 which removes the effect of the interleaver 206 at the 
transmitter. 

5 

The decoder 230 maps iteratively decodes the stream of bit level LLR's as 
a soft output by use of a soft input soft output (SISO) algorithm, such 
algorithms include Maximum-A-Posteriori (MAP) or Low Density Parity 
Check (LDPC) algorithms. 

10 

Optionally, information obtained from the decoder 230 at each iteration is 
fed back into the conditional probability calculations, typically this is 
symbol probability information. This information is treated as a-priori 
information which increases the accuracy of the information obtained 
15 according to: 



P(Cj =S p )*_ £ fl exp 

C:Cj=S p i=x+\-E c 

V J 

After a number of iterations, as few as ten have been shown to yield 
improvements in performance over prior art systems, final, hard, 
definitive, bit decisions are based upon the soft output values from the 
20 decoder 230. 

The method detailed hereinbefore estimates which symbol is most likely 
to have been transmitted using the Maximum-A-Posteriori (MAP) criteria 
in which: 



Zi R i,k C k 



m 



/=*+! 



2\ 



2a 



m 



np r [c s ] 

s=x+\-Ec,s*j 
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A - 

Cj = arg max V(C:Cj=S p ) Pr (C : Cj = S p ) 

S f 6[C t C, C^J 

Referring now to Figure 5, a method of signal processing for a MIMO 
system comprises defining a QR decomposition of a channel gain matrix 
and setting an enumeration constraint upon the matrix (Step 500). 
5 Conditional symbol probabilities are calculated by enumerating across the 
desired block of the R matrix, using the enumeration constraint (Step 
502). The effect of the symbol having the highest conditional symbol 
probability within the desired block of the R matrix is cancelled out (Step 
504) and the R matrix is traversed with the cancellation process being 
10 repeated for each enumeration block (Step 506). 

9 

The apparatus and method described hereinbefore reduce the 
computational complexity of symbol determination by trading diversity 
upon receive against computational complexity. 

15 

Using the MAP criterion to form decisions as in this approach gives a 
retained diversity order in symbol estimation that is given by the 
enumeration constraint E c . 

20 Cj = arg max «F (C:Cj=S p ) Pr (C : Cj = S p ) < MAP Criterion) 

V<c,c 2 c 2mc ) 

The constraint E c establishes smaller areas for symbol enumeration 
instead of using the whole channel matrix. Larger areas lead to a larger 
list of possible symbol candidates. Computational complexity is 
25 maintained within areas of equal size as the matrix is traversed. 
Enumeration are still related to the modulation mode utilised and is given 

* 

by: 
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The approach provides a flexible LLR computation, which is not 
restricted to the number of transmit antennas used in the MIMO 
5 architecture but upon a predefined subset of the transmit antennas at the 
receiver. The use of feedback within the computational scheme results in 
an improved computational scheme in which detection and decoding. are 
not disjointed. This is of importance as any improvement in detection 
yields an improvement in decoding results and vice-versa. 
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CLAIMS 

1. A signal processing apparatus comprising a plurality of receiving 
5 means arranged to receive a composite signal indicative of a plurality of 

symbols transmitted, substantially simultaneously, from a plurality of 
remote transmission means, and processing means arranged to iteratively 
decode each most probable symbol contained in said composite signal, 
within a constrained enumeration formalism. 

10 

2. Apparatus according to Claim 1 wherein the processing means is 
arranged to define an enumeration constraint for use in the constrained 
enumeration formalism. 

15 3. Apparatus according to either of Claims 1 or 2 wherein the 
processing means is arranged to perform a QR decomposition upon a 
channel gain matrix. 

4. Apparatus according to Claim 3 wherein the enumeration constraint 
20 is a number of entries in an R matrix over which probable symbols are 

* 

enumerated. 

5. Apparatus according to any preceding claim wherein the processing 
means is arranged to determine the most probable symbol by enumerating 

25 across all symbol conditional probabilities for each possible symbol. 

6. Apparatus according to Claim 5 wherein the processing means is 
arranged to convert symbol conditional probability to a bit level 
logarithmic likelihood ratio (LLR). 

30 
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7. Apparatus according to any preceding claim including a parallel to 
serial conversion means arranged to convert parallel, bit level, LLR's 
into a single stream of LLR's. 

5 8. Apparatus according to Claim 7 including a deinterleaving means 
arranged to deinterleave, bit level, LLR's from the single stream of 
LLR's. 

■ 

9. Apparatus according to any preceding claim including decoding 
10 means arranged to apply iterative soft input soft output (SISO) decoding 

to single bit LLR's to determine a symbol. 

10. Apparatus according to Claim 9 wherein the decoding means is 
arranged to pass a symbol probability to the processor for inclusion in an 

15 iterative enumeration step. 

11. Apparatus according to either of Claims 9 or 10 including a hard 
decision unit that is arranged to determine a symbol based upon a soft 
output from the decoding means. 

20 

12. A method of signal processing for a MIMO system comprising the 
steps of: 

i) receiving a composite signal indicative of a plurality of symbols; 

♦ 

ii) performing a QR decomposition upon a channel gain matrix for 
25 the composite signal; 

iii) defining an enumeration constraint; 

iv) calculating possible conditional probabilities for one of the 
plurality of symbols contained within the composite signal, using 
the enumeration constraint; and 
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v) iterating step iv), incorporating a most probable symbol for the one 
symbol determined in the previous iteration of step iv) in the 
conditional probability calculation operation. 

5 13. A method according to Claim 12 including setting the enumeration 
constraint to encompass a sub-set of possible transmit antennas. 

14. A method according to either of Claims 12 or 13 including defining 
the enumeration constraint as a number of elements within an R matrix. 

10 

15. A method according to any one of Claims 12 to 14 including 
calculating a symbol conditional probability in order to determine the 
most probable symbol received over a given transmission channel. 

15 16. A method according to Claim 15 including converting the symbol 
conditional probability to a bit level logarithm likelihood ratio (LLR). 

17. A method according to any one of Claims 12 to 16 including 
converting a plurality of parallel streams of bit level LLR's to a serial 

20 stream of bit level LLR's. 

18. A method according to Claim 17 including deinterleaving bit level 
LLR's from the serial stream of bit level LLR's. 

25 19. A method according to any one of Claims 16 to 18 including 
decoding the single bit LLR's. 

20. A method according to any one of Claims 12 to 19 including 
making a hard determination of a received symbol based upon a soft 
30 output from the decoding operation. 
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21. A method of reducing the computational load of a signal processor 
in MIMO architectures comprising the steps of: 

« 

i) receiving composite input signals having spatial diversity from each 
of a set of n receiver elements; 
5 . ii) constructing an n by m channel matrix from values indicative of 

channel gains between each transmit and receive element; 

iii) executing a QR decomposition upon the channel matrix to form an 
upper triangular R matrix and a unitary Q matrix; 

iv) enumerating to determine probabilities of a given symbol being 
transmitted from a given transmitter using a constrained data sub-set 
of the triangular matrix; and 

v) making a hard decision about which possible symbol is the most 
probable symbol to have been transmitted so as to reduce the number 
of enumerations required to carry out a further probability 
calculation. 

22. The method of Claim 21 including using sub-optimally determined 
symbol values to generate final definite symbol values. 

23. A computer readable medium having therein instructions for causing 
a processing unit to execute the method of any one of Claims 12 to 20. 

24. A computer readable medium having therein instructions, for causing 
a processing unit to execute the method of either of Claims 21 or 22. 
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